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Preface 

We  live  in  a  very  complex  world.  Time  was  when  a  man 
could  be  a  jack-of-all-trades  and  master  them  all.  Not  so 
anymore.  Today  a  jack-of-all-trades  would  be  master  of  none. 
Becoming  a  master  of  just  one  trade  now  requires  near  total 
committment  to  it  not  just  in  the  apprenticeship  stages,  but 
also  in  the  years  that  follow  merely  to  keep  abreast  of  the 
constant  changes. 

Having  been  a  full-time  U.  S.  Air  Force  pilot  for  the 
six  years  prior  to  arriving  at  the  Air  Force  Institute  of 
Technology,  I  found  out  the  truth  of  that  statement  very 
quickly.  I  earned  my  B.  S.  degree  in  electrical  engineering 
six  years  before,  but  by  this  time  neglect  of  all  things 
mathematical  had  taken  its  toll  on  any  engineering  skill 
I  may  once  have  possessed.  So,  my  19-month  tenure  at  AFIT 
was  a  long,  hard  struggle  to  re-establish  myself  in  the 
engineering  world. 

It  is  for  this  reason  that  I  am  truly  grateful  to  my 
thesis  advisor.  Major  James  D.  Dillow,  for  his  constant 
encouragement  and  help  on  the  work  presented  herein.  I 
would  also  like  to  thank  Mr.  Ronald  0.  Anderson  of  the 
Flight  Simulation  Branch,  Air  Force  Flight  Dynamics  laboratory, 
for  his  support  on  this  project.  A  special  thanks,  also, 
to  Mr.  Allan  Carney  for  his  eager  assistance  in  setting  up 
the  analog  simulation.  Upon  his  aid  much  of  the  success 
of  the  simulation  rests. 
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Of  course,  this  work  would  never  have  gotten  off  the 
ground  without  the  cooperation  of  my  three  pilots.  I  am 
deeply  indebted  to  Captains  Jim  Butt,  Ed  Jordan,  and  Jim 
Tilley  for  foregoing  what  precious  little  free  time  students 
at  AFIT  enjoy  and  spending  many  hours  apiece  flying  my 
simulation. 

And,  finally,  for  putting  up  with  a  non-existent 
husband  for  these  past  19  months,  a  special  debt  of 
gratitude  is  extended  to  my  wife,  Trish.  She  bravely 
carried  on  through  thick  and  thin. 


Thomas  R.  Harvey 
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Abstract 

The  two-dimensional  (longitudinal)  kinematics  of 
air-to-air  combat  tracking  with  a  lead  computing  optical 
sight  system  were  simulated  on  an  analog  computer.  Three 
pilots  flew  three  different  aircraft  configurations  on  the 
fixed-base  simulator  at  3000-  and  1000-foot  ranges  against 
a  target  driven  to  RMS  accelerations  of  3.5  G  and  5.0  G 
by  filtered,  white,  Gaussian  noise.  Averaged  RMS  data 
were  recorded  on  attacker's  elevator  deflection,  pitch 
rate,  lead  angle,  line  of  sight  to  the  target,  and  tracking 
error  for  each  case. 

Simultaneously,  the  identical  tasks  were  performed  by 
a  digital  pilot  model  formulated  from  optimal  control  theory 
with  visual  acuity  threshold  effects  accounted  for.  A 
comparison  of  the  data  generated  by  the  human  pilots  versus 
that  of  the  pilot  model  showed  very  nearly  a  one-to-one 
correlation  for  elevator  deflection,  lead  angle,  and  line 
of  sight.  However,  the  results  were  not  quite  as  comparable 
when  applied  to  pitch  rate  and  tracking  error.  A  re- 
evaluation  of  cost  functional  weightings  based  upon  a 
further  analysis  of  pilot  response  in  this  task  led  to 
model  performance  that  correlated  with  the  human  pilot 
performance  data  very  well. 
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APPLICATION  OF  AN  OPTIMAL  CONTROL 
PILOT  MODEL  TO  AIR-TO-AIR  COMBAT 


I.  Introduction 


A  Problem  Exists 

All  too  often  production  engineering  runs  into  design 
problems  at  the  final  stages  of  production  due  to  incom¬ 
patibility  of  integrated  subsystems.  This  problem  is 
especially  critical  when  the  system  is  being  designed 
for  use  by  a  human  operator.  The  reason  for  this  anomaly 
is  the  fact  that,  even  though  individual  subsystems  may  be 
designed  perfectly  for  the  task  that  is  to  be  performed, 
when  they  are  integrated  by  means  of  a  human  operator,  the 
overall  system  can  break  down  because  of  the  operator's 
inherent  complexities  and  imperfections.  For  example,  if 
a  subsystem  is  designed  with  very  little  stability  margin, 
the  lag  time  inherent  in  the  human  operator's  reaction  might 
well  cause  the  system  to  become  unstable.  This  is  due  to 
the  fact  that  the  action  of  the  operator  could  actually  de¬ 
grade  the  combined  action  rather  than  achieve  the  stabilizing 
effect  intended.  Therefore,  increased  emphasis  is  being 
placed  on  total  system  simulation  with  the  pilot  in  the  loop 
during  the  preliminary  design  stages  of  production. 

These  system  simulations  usually  take  one  of  two 
forms.  The  first  is  direct  simulation  of  the  total  system 
on  analog  or  hybrid  computers.  This  demands  a  human  oper¬ 
ator  to  close  the  control  loop  by  performing  the  design  task. 
The  second  method  involves  mathematical  simulation  of  the 
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operator.  This  model  is  used  to  predict  human  operator 

t 

reaction  to  system  design  and  to  evaluate  the  compatibility 
of  the  integrated  system.  The  latter  method  alleviates  the 
need  for  the  repetitive  data  collection  process  required  when 
the  human  operator  is  used  in  the  simulation,  so  long  as  it 
can  be  shown  that  the  mathematical  pilot  model  successfully 
duplicates  human  pilot  reaction.  The  objective  of  this 
analysis  is  to  show  that  a  mathematical  pilot  model  can 
accurately  predict  the  performance  of  integrated  aircraft 
flight  control  and  lead  computing  optical  sight  (LCOS) 
systems  in  performing  the  air-to-air  tracking  task,  fighter- 
on-  fighter. 

Most  studies  in  the  past,  such  as  the  many  "Paper  Pilot" 

(Ref  1)  investigations  have  utilized  the  quasi-linear  pilot 
describing  function  models  as  formulated  by  McRuer  (Ref  14) . 

However,  the  present  study  uses  the  optimal  pilot  describing 
function  model  formulated  by  Kleinman,  Baron,  Miller,  Elkind, 
and  Levison,  hereafter  referred  to  as  the  "Kleinman  model" 

(Ref  2) .  This  model  incorporates  a  Kalman  filter  and  state 
predictor  to  simulate  the  analytical  and  anticipatory  tasks 
which  a  pilot  must  perform  to  properly  control  an  aircraft. 

Kleinman  has  shown  th.t  this  model  very  accurately  represents 
the  response  of  a  pilot  in  performance  of  a  number  of  tasks, 
such  as  the  landing  approach  (Ref  8) ,  helicopter  hover  (Ref  2) , 
and  ground-based  anti-aircraft  gunnery  (Ref  10) . 

The  Analysis  Procedure 

To  show  that  the  Kleinman  pilot  model  does  indeed  exhibit 
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performance  characteristics  closely  matched  to  those  of  a 
human  pilot  in  the  air-to-air  tracking  task,  two  separate 
simulations  are  required.  First,  data  must  be  gathered  on 
actual  pilot  performance  in  this  task.  That  is  accomplished 
in  this  study  by  simulating  the  task  on  an  analog  computer 
and  presenting  it  to  three  pilots  with  previous  fighter  air¬ 
craft  experience.  This  presentation  is  done  via  a  fixed-base 
simulator  to  be  described  in  Section  III.  The  data  thus 
collected  reflect  pilot  workload  and  performance  of  the  task 
confronting  him  in  each  experimental  case. 

The  second  simulation  pits  the  Kleinman  model  against 
the  identical  tasks  on  the  CDC  6600  digital  computer.  The 
same  statistical  data  on  workload  and  performance  are  collected 
from  the  simulation  and  compared  with  that  obtained  from  the 
analog  simulation.  The  results  of  this  comparison  are  re¬ 
vealed  in  Section  V  of  this  report. 

To  obtain  a  large  enough  spread  in  the  data  to  accentuate 
the  correlation  between  analog  and  digital  results,  several 
different  aircraft,  ranging  from  good  to  bad  aerodynamic 
designs  for  this  task,  are  used  at  varying  altitudes  and 
Mach  number.  Additionally,  for  each  aircraft,  two  different 
tracking  ranges  and  two  different  target  RMS  ”G"  loadings 
are  simulated. 

Limits  and  Assumptions 

The  dynamics  of  the  air-to-air  combat  tracking  task  are 
extremely  complex  when  considered  from  a  three-dimensional 
aspect.  This  is  due  primarily  to  the  inertial  cross-coupling 
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of  the  aircraft  longitudinal,  lateral,  and  yaw  axes  that 
manifests  itself  in  the  aircraft  equations  of  motion.  So, 
as  a  starting  point  in  the  analysis  of  this  problem,  this 
study  is  limited  to  an  analysis  of  the  longitudinal  tracking 

task  only,  and  neglects  any  cross-coupling  effects.  This 

\ 

procedure  is  common  practice  and  usually  lends  itself  to 

\ 

realistic  results.  Future  work  should  concentrate  on 
expansion  of  this  study  to  three  dimensions. 

It  is  further  assumed  that  the  problem  begins  with  the 
attacking  aircraft  already  in  a  tracking  position  relative 
to  the  target  aircraft;  thus  the  task  is  "merely"  to  continue 
tracking  the  target  as  closely  as  possible.  The  attacking 
aircraft  is  considered  to  be  rigid  and  impervious  to  "G" 
loadings  in  excess  of  normal  maximums  in  the  interest  of 
tracking  the  target  with  minimum  error.  Also,  it  is 
assumed  that  the  attacking  aircraft  possesses  the  capability 
of  producing  infinite  thrust  and  can  therefore  maintain 
constant  Mach  number  under  any  conditions  of  G- loading, 
angle  of  attack,  or  attitude.  Since  the  Kleinman  pilot 
model  is  basically  a  linear  regulator,  all  dynamic  equations 
in  the  study  have  been  linearized  about  the  aircraft  sta¬ 
bility  axis,  and  perturbations  are  assumed  small  enough  to 
validate  the  use  of  these  linear  equations  in  a  nonlinear 
environment.  In  addition,  the  aircraft  stability  derivatives 
are  assumed  constant  throughout  the  tracking  task,  an  assump¬ 
tion  which,  in  the  light  of  assumed  constant  Mach  number,  is 
quite  valid. 

i 

Finally,  the  affects  of  gravity  on  this  problem  contribute 
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only  bias  terms  that  merely  complicate  the  equations  and 
computations .  Therefore ,  these  affects  are  neglected, 
and  any  reference  to  "6”  forces  in  this  text  alludes  only 
to  kinematic  accelerations. 
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II.  The  Air-To-Air  Tracking  Task 

Modelling  the  Sight  Dynamics 

The  task  of  accurately  tracking  a  maneuvering  target 
for  the  purpose  of  obtaining  a  "kill"  with  airborne  cannon 
fire  is  one  of  the  most  difficult  required  of  a  fighter 
pilot.  As  is  the  case  when  firing  any  projectile  at  a 
moving  target,  lead  for  target  motion  must  be  computed 
and  the  aiming  direction  adjusted  to  compensate  accordingly. 
Further  complicating  the  air-to-air  problem  are  trajectory 
adjustments  that  must  be  made  to  account  for  projectile 
drag,  velocity  jump,  and  gravity  drop,  as  well  as  other 
more  minor  effects.  The  total  compensation  is  normally 
presented  to  the  pilot  of  the  attacking  aircraft  in  the 
form  of  a  two-mils-in-diameter  "pipper",  surrounded  by 
a  larger  circle  called  the  reticle,  and  depressed  from  the 
weapon  line  by  the  amount  of  this  computed  lead  angle,  as 
is  shown  in  Figure  1.  Thus,  when  the  pilot  maneuvers 
his  aircraft  in  such  a  manner  as  to  place  the  "pipper" 
on  the  target,  he  has  achieved  the  proper  aiming  direction 
to  insure  a  kill  should  the  target  maneuver  remain  constant. 

Neglecting  windage  jump  and  the  effects  of  gravity, 
the  lead  angle  required  is,  in  general,  computed  from  the 
following  vector  equation  (Ref  18:7-8): 

T  V 

X  =  tT  +  (s  x  VT)  -  Jy  o  (1) 

where 


X  »  vector  lead  angle 
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D  *  present  range 

V,  =  mean  relative  velocity  of  projectile  over  one 
time  of  flight 

=  inertial  angular  rate  of  line  of  sight  (LOS) 
from  the  attacker  to  the  target 

=  time  of  flight  to  target  future  position 

s  =  unit  vector  along  the  LOS 

V,^  =  vector  acceleration  of  target  along  vertical 
(positive  down) 

Jv  *  ballistic  jump  parameter 
VA  =  attacke-;  velocity  magnitude 

a  =  angle  of  attack 

The  first  term  in  this  equation  is  that  portion  of  the 
lead  angle  that  compensates  for  linear  target  motion.  The 
second  term  accounts  for  target  acceleration  across  the 
line  of  sight,  the  vector  cross-product  being  the  component 
of  target  acceleration  normal  to  the  LOS.  The  last  term 
adjusts  the  lead  angle  for  initial  trajectory  errors  due 
to  attacker  angle  of  attack  at  the  instant  of  firing. 

Figure  2  portrays  the  principle  angle  relationships  in 
the  in-plane,  air-to-air  tracking  situation.  Since  the  LOS 
rate  (£T)  cannot  be  directly  measured  by  the  attacking  air¬ 
craft,  it  is  approximated  in  the  LCOS  system  by  ,  the  . 

inertial  angular  rate  of  the  reticle,  or  pipper.  This 
approximation  is  good  so  lonj  as  the  sight  is  held  exactly 
on  the  target.  However,  when  the  attacker  is  maneuvering 

to  put  the  sight  on  the  target,  or  the  target  does  not 

•  • 

remain  in  a  constant  maneuver,  Z  is  not  equal  to  I_,  and 

r  ^ 

the  computed  lead  angle  will  be  in  error  until  steady 
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This  is  the  equation  that  is  instrumental  in  the  LCOS 
system  to  position  the  pipper  at  the  current  correct  aiming 
point  on  the  pilot's  "heads-up"  display. 


Velocity  Jump  Correction 

Jy*  th'*  velocity  jump  correction  angle,  is  illustrated 
in  Figure  3.  It  is  a  consequence  of  the  fact  that  the 
projectile  muzzle  velocity  vector  and  the  aircraft  velocity 
vector  are  not  normally  colinear  when  the  aircraft  is  being 
maneuvered,  and  is  calculated  from  (Ref  20:47) 


^A  x  —  (WL) 


V.  sin  a 


where  I(«nj 


unit  vector  directed  along  the  weapon  line. 


Assuming  small  angles 


Jv  S  V  +~V  (22) 

A  m 

Because  it  is  a  function  of  angle  of  attack  (a) , 
the  jump  correction  function  had  to  be  linearized  about  a 
nominal  angle  of  attack  in  order  to  render  equation  (20) 
linear.  This  nominal  value  was  chosen  as  10°  in  all  cases. 


Short-Comings  of  the  LCOS  System 

Once  the  attacker  is  in  a  tracking  position  on  the 
target  (normally  behind  and  slightly  elevated  in  the  target's 
plane  of  motion) ,  he  must  then  generate  an  angular  rate  of 
turn  in  that  plane  approximately  equal  to  that  of  the  target. 
The  major  tracking  problem  with  the  LCOS  system  arises  when 
the  attacker  does  not  have  the  sight  on  the  target.  For 
instance,  if  the  sight  is  behind  the  target  the  pilot  must 
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"V 


increase  his  turn  rate  to  catch  up.  This  results  in  an 
erroneous  increase  in  the  displayed  lead  angle,  leaving 
the  pilot  with  no  valid  information  as  to  exactly  where  to 
position  his  aircraft  until  his  rate  of  turn  is  once  again 
constant  and  equal  to  that  of  t  e  target.  A  process  of 
"hunting"  then  ensues  until  the  pilot  has  either  found  the 
exact  solution  or  decides  to  hold  one  solution  steady  and 
fire  continuously  while  allowing  the  sight  to  "slide  through" 
the  target.  In  practice,  the  latter  procedure  usually  leads 
to  the  best  results.  As  is  evident  from  equation  (20) ,  as 
range  increases,  the  sight  time  constant,  T^,  becomes  larger, 
causing  these  unwanted  dynamics  to  become  even  worse. 

Normally,  the  sight  is  designed  with  a  nominal  time  constant 
based  upon  a  nominal  firing  range  of  around  1500  feet  (Ref  5:1-1). 

i 


13 


GA/MA/74M-1 


V 


In  Appendix  A  it  is  pointed  out  that  the  sight  tine  con¬ 
stant  for  each  case  .in  this  study  was  computed  based  upon 
the  given  tracking  conditions. 


Kinematics 

The  inertial  line  of  sight  from  the  attacker  to  the 
target  is  computed  from  (Ref  18:11) 

lr*5  lS*vT-£*v*l  <23> 


which  states  that  the  inertial  rate  of  rotation  of  the  LOS 
is  given  by  the  difference  between  the  relative  angular 
rates  of  rotation  of  the  target  and  the  attacker.  From 
Figure  2  this  equation  becomes,  considering  the  in-plane 
task  only 

V  V 

zT  ■  ir  sin(*T  **  V  “  it  sin<x-a-e>  (24> 

Using  small  angle  assumptions  and  noting  from  Figure  2 

that 

e  =  IT  -  Ep  =  IT  -  0-  +  X  (25) 

V  V 

=  {lfT  ’  V  "  (6"a‘ZT>  <26) 


or 


,  _  VA 

£T  D-  a 


V  -V  V  V 

+  A  VT  .  _T  Y  _  _A 

+  D  £T  +  D  YT  D  0 


(27) 


Since  allowing  a  rate  of  change  of  range,  D,  would  render 
the  equations  nonlinear,  the  cases  considered  in  this  study 
were  limited  to  zero  closure  rate  with  the  velocities  of 
the  attacker  and  target  constant  and  equal.  Therefore,  the 
second  term  in  (27)  was  always  zero  in  this  study. 
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III.  Analog  Simulation  of  the  Air-to-Air  Task 

The  purpose  of  the  analog  computer  simulation  was  to 
gather  appropriate  data  on  specific  performance  functions 
to  compare  with  similar  results  from  digital  computer  runs 
of  the  identical  tasks  "flown"  by  the  analytical  pilot 
model.  A  system  model  of  this  task  can  be  broken  up  into 
five  functional  areas,  as  shown  in  Figure  4:  the  pilot  (or 
pilot  model) ,  the  flight  control  response,  the  aircraft 
equations  of  motion,  the  sight  computer,  and  the  kinematics 
of  the  problem.  The  system  variables  chosen  as  the  states 
which  fully  describe  the  task  were  the  attacking  aircraft's 
elevator  deflection,  6,  pitch  angle,  6,  pitch  rate,  q,  and 
angle  of  attack,  a;  the  lead  angle,  X,  the  inertial  line  of 
sight  from  the  attacker  to  the  target,  ZT;  the  normal 
acceleration  of  the  target,  A  ;  the  flight  path  angle 
of  the  target,  y T?  and  a  dummy  state  representing  the 
first  stage  filtering  of  white,  Gaussian  noise. 

The  Attacker  Model 

Equations  of  Motion.  As  was  mentioned  previously  in 
the  introduction  to  this  text,  the  equations  of  motion  of 
the  attacking  aircraft  were  approximated  by  the  standard 
"short-period"  stability  axis  equations.  These  equations 
are  first-order,  linear,  and,  in  view  of  the  assumption 
of  constant  stability  derivatives,  time-invariant  differential 
equations  which  approximate  aircraft  response  to  control  inputs. 
Since  this  analysis  is  concerned  only  with  the  longitudinal 
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Figure  4,  Functional  Diagram  of  the  Two-Dimensional  Air-To-Air  Tracking  Task. 


GA/MA/74M-1 


axis  dynamics,  only  three  equations  are  necessary.  Those 
describing  pitch  angle  and  angle  of  attack  were  presented 
in  Section  II  as  equations  (10)  and  (11),  respectively. 
The  third  equation,  that  for  pitch  rate  is  (Ref  4:14-21) 

q**Mq  +  Ma  +  M*a  +  M.6  (2( 

q  a  c  o 

substituting  (11)  in  (28)  and  rearranging 

q  =  (M.+M-Z  .)£  +  (M  +M-)q  +  (M  +M*Z  )o  (2 S 

6  a  6  q  a  a  a  a 

Elevator  Actuation.  The  transfer  function  relating 

elevator  deflection,  £,  to  commanded  elevator  deflection, 

£  ,  for  a  typical  tactical  fighter  aircraft  is  given  by 

(Ref  19:155) 


£  1  +  T  s 

c  a 

where  KT  is  the  control  linkeage  gain  and  x  is  the  time 
constant  of  the  power  cylinder  actuator.  Typical  values 
for  these  two  parameters  are  0.8  (dimensionless)  and  0.05 
seconds,  respectively.  The  differential  equation  for 

elevator  deflection  is  then 

1  Kt 

$  =  _  £  +  _id  £c  (31) 

a  a 

The  analog  simulation  utilizes  a  force  stick  that 
converts  pilot  pressure  in  pounds  to  volts,  where  one  volt 
equals  0.01  radian  in  the  scaled  circuitry.  The  equation 
modelled  in  the  circuitry  is  (Ref  3:12) 


-  £  =  K,F 

c  f  s 


where  Kf  is  the  force  stick  sensitivity  in  units  of  radians 
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per  pound,  and  F  is  the  applied  force  in  pounds.  The 

S 

negative  sign  on  6  is  a  consequence  of  the  fact  that 

C 

positive  elevator  deflection  is  commonly  defined  in  the 
aircraft  equations  of  motion  as  that  movement  which  results 
in  a  negative  pitch  rate;  i.e.,  "down”  elevator.  Sub¬ 
stitution  of  (32)  into  (31)  then  gives 

j  .  .  i_  4  +  ,  ,33, 

a  a 

Figure  6  is  the  analdg  schematic  of  the  attacking 
aircraft  dynamics.  Note  that  a  -stick  sensitivity  adjustment 
potentiometer  was  included  in  the  factors  making  up  K^. 

This  "pot"  was  located  on  the  instrument  panel  of  the 
simulated  aircraft  cockpit  so  that  the  pilot  could  adjust 
the  stick  sensitivity  to  an  optimum  level  for  his  technique. 
This  was  done  so  that  limitations  on  stick  "feel"  would  not 
be  a  factor  in  each  pilot's  performance  of  the  task. 

Target  Modelling 

Generation  of  Target  Motion.  The  target  aircraft  was 

a  fighter  aircraft  possessing  capabilities  at  least  as 

strong  as  those  of  the  pursuer.  Target  normal  acceleration, 

A  ,  was  derived  from  band-limited,  Gaussran  white  noise 
nT 

with  statistics  of  zero  mean  and  standard  deviations 
(RMS)  of  3.5G  ar.d  5.0G  for  each  case  tested.  The  analog 
schematic  of  Figure  6  was  used  to  record  target  motion 
for  purposes  of  repetitive  use.  The  output  of  a  Gaussian 
white  noise  generator  was  subjected  to  a  double  filter 
network  to  obtain  a  random  function  which  represented 
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Figure  6.  Analog  Circuit  for  Taping  Target  Motion 


"realistic"  target  motion.  The  break  frequency,  1/tt, 
was  manipulated  by  trial  and  error  until  this  "realistic" 
motion  was  obtained.  A  value  of  tt  equal  to  3  seconds 
was  found  to  be  best  for  this  situation.  It  was  quickly 
discovered  that,  even  though  the  tape  recorded  "motion"  of 
the  target  was  completely  random,  characteristic  peaks  in 
the  function  were  easily  remembered  by  the  test  pilots,  and, 
after  a  few  runs,  even  anticipated.  To  alleviate  this  prob¬ 
lem,  three  separate  target  motions  were  recorded  and  used 
in  each  case.  Figure  7  shows  reproductions  of  the  three 
target  motions  used. 

Target  Interface  with  the  System.  Recall  from  Section 
II  that  equation  (27)  called  for  a  knowledge  of  the  target 
flight  path  angle,  It  should  be  noted  here  that  this 

parameter  is  used  explicitly  in  this  study  only  because  it 
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is  available  in  the  simulation.  In  reality  it  would  be 
approximated  by  own  ship's  flight  path  angle.  However, 
since  it  is  available  in  the  simulation  and  since,  as  with 
equation  (7) 


then 


yTVT 


Yt 


(34) 


(35) 


Figure  8  is  the  analog  schematic  of  target  motion 
input  to  the  system  dynamics,  specifically,  the  £T  equation. 
The  bias  circuitry  was  necessitated  by  the  requirement  for 
zero  mean  motion  (see  Section  IV) ,  and  the  gain  adjustment 
was  to  facilitate  changing  RMS  "G"  levels.  The  "1/VT"  term 
in  (35)  forced  normalization  of  the  equations  since  that  term 


Figure  8.  Analog  Circuit  for  Target  Flight 
Path  Angle  Generation 
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the  target  signal  from  reaching  the  rest  of  the  system.  One 
thousand  feet  was  chosen  as  one  distance  unit  (DU)  in  the 
normalized  equations,  with  one  second  remaining  one  time 
unit  (TU) .  Consequently,  one  velocity  unit  (VU)  became 
1000  ft/sec  and  one  acceleration  unit  was  1000  ft/sec2. 
Though  it  does  not  appear  in  the  equation,  a  dummy  feedback 
path  was  added  to  the  yt  equation  to  keep  the  steady  state 
covariance  of  yt  from  becoming  infinite  (see  Section  IV) . 

A  nominal  value  of  0.01  was  used  for  this  feedback  path 
so  that  it  would  have  little,  if- any,  affect  on  the  actual 
dynamics  beind  modelled. 

Sight  Dynamics  and  Problem  Kinematics 

The  lead  angle  and  line  of  sight  angle  equations 

were  derived  in  Section  II  as  equations  (20)  and  (27) , 

respectively.  The  analog  models  of  these  equations  are 

shown  in  Figure  9.  Also  shown  are  the  circuits  for  the 

computation  of  I  ,  the  relative  line  of  sight  of  the 

A 

target  with  respect  to  the  attacker  weapon  line,  and  e, 
the  tracking  error.  The  mathematical  relationships  for 
these  last  two  parameters  were  derived  from  the  geometry 
of  Figure  2  as 


and 

e  »  X  -  ZT  (37) 
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Figure  9.  Analog  Representation  of  Sight  and 
Problem  Kinematics 

The  Fixed-Base  Simulator 

Figure  10  is  a  photograph  of  the  simulator  "cockpit". 
It  consists  of  a  fighter  aircraft  ejection  seat  with  side- 
stick  controller  and  a  dual-beam  cathode  ray  oscilloscope 
upon  which  the  task  is  displayed  mounted  at  eye  level.  The 
side-stick  controller  is  similar  to  that  being  used  in  one 
advanced,  light-weight  fighter  "fly-by-wire"  design.  In 
the  simulation,  the  stick  is  trimmed  for  level  flight  and 
cannot  be  retrimmed. 

Presentation  of  the  Task  to  the  Pilot.  When  tracking 
a  target,  the  pilot  of  the  pursuing  aircraft  is  aware  prima 
rily  of  two  parameters:  the  position  of  the  target  and  the 
position  of  the  sight.  The  difference  between  the  two 
positions  is  the  tracking  error,  and  this  is  what  the  pilot 


attempts  to  minimize.  In  this  simulation,  the  target  is 


GA/MA/74M-1 


displayed  on  the  oscilloscope  as  an  inverted  "T".  It  is 

positioned  relative  to  the  center  of  the  oscilloscope, 

assumed  to  be  the  extension  of  the  aircraft's  weapon  line, 

by  the  current  magnitude  of  E_  ,  the  relative  line  of  sight 

A 

to  the  target.  The  sight  is  positioned  relative  to  the 
center  of  the  oscilloscope  by  the  current  magnitude  of  X, 
the  lead  angle.  Figure  11  illustrates  the  picture  seen  by 
the  pilot.  The  tracking  error  shown  is  positive;  i.e.,  the 
attacker  is  lagging  behind  the  target.  Note  that  only  the 
reticle  portion  of  the  sight  is  represented.  This  is  due 
to  the  limitations  encountered  on  being  able  to  physically 
generate  both  the  reticle  and  the  pipper,  as  well  as  the 
target  symbol,  on  the  same  oscilloscope  at  the  same  time. 
The  pilots  thus  had  to  estimate  the  position  of  the  center 
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act  ended 
weapon  llne'X 


reticle 


Figure  11.  Simulation  Display 


of  the  reticle  to  achieve  zero-error  tracking. 

The  Gathering  of  Data.  After  an  extended  period  of 
training  during  which  the  three  pilots  became  used  to  the 
task  before  them  and  the  feel  of  the  side-mounted,  force- 
stick  controller,  the  taking  of  data  began.  Twelve  separate 
cases  were  run  using  three  different  aircraft  dynamics,  two 
ranges,  and  target  RMS  G  levels  of  3.5G's  and  5.0G's.  The 
three  aircraft  were  the  F-4E,  the  F-5,  and  the  A-7.  A 
summary  of  the  dynamics  of  each  aircraft  is  presented  in 
Table  I.  Ranges  of  3000  feet  and  1000  feet  were  decided  upon 
since  they  represent  approximate  maximum  and  minimum  tracking 
ranges  in  an  actual  situation.  As  was  mentioned  previously, 
three  different  recorded  target  motions  possessing  identical 
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Table  I. 

Dynamic  Parameters  and  Stability  Derivatives 


Aircraft 

F-4Et 

F-5* 

|  A-7* 

Altitude  (Ft) 

15,000 

5,135 

15,000 

Mach  Number 

0.90 

0.81 

0.60 

Velocity,  VQ  (Ft/Sec) 

951.6 

889.0 

635.0 

Dynamic  Pressure,  qQ 

(Lb/Ft2)  677.3 

804.4 

301.0 

Mass,  m  (Slugs) 

1433.5 

354.0 

680.0 

Reference  Area,  S  (Ft 

2)  530.0 

170.0 

375.0 

C  (1/Rad) 

— 

3.95 

4.40 

a 

CT  (1/Rad) 

L6 

-  - 

0.859 

0.600 

Z  **  (1/Sec) 
a 

-1.0326 

-1.7164 

-1.1502 

Z6**  (1/Sec) 

-0.09512 

-0.3733 

-0.1568 

Mo  (1/Sec2) 

-10.443 

-10.30 

-9.08 

(1/Sec) 

-0.3439 

-0.0646 

-0.133 

M  (1/Sec) 

-0.7381 

-1.350 

-0.696 

(1/Sec2) 

-37.08 

-47.2 

-18.90 

t  (Ref  7s 80-81) 

*  (Ref  17:146,148) 

**  Values  for  Z  and  Z,  for  the  F-5  and  the  A-7  were 
a  o 

calculated  from 
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statistics  were  used  in  each  case.  Thus,  each  pilot  "flew” 
a  minimum  of  36  times  once  data-gathering  runs  began. 

A  typical  data  run  began  with  the  target  and  the  sight 
both  stationary  in  the  center  of  the  screen,  making  all 
initial  conditions  equal  to  zero.  The  picture  thus  dis¬ 
played  represented  a  situation  in  which,  unbeknown  to 
the  target,  the  attacker  had  maneuvered  into  a  stern,  or 
"6  o'clock”,  attack  position.  When  the  pilot  was  ready, 
a  switch  was  thrown  to  activate  the  problem  dynamics  with 
no  target  motion  other  than  that  generated  by  the  action 
of  the  attacker  pilot.  Shortly  thereafter,  a  second  switch 
was  thrown  sending  the  taped  target  motion  to  the  display. 
This  simulated  sudden  target  awareness  of  the  attacker's 
approach  followed  by  high  "G"  evasive  action.  At  a  pre¬ 
determined  point  on  the  tape  a  switch  was  thrown  to  start 
the  data  measurement  circuitry.  For  100  seconds  thereafter, 
analog  integrators  evaluated  the  time  integrals  of  elevator 
deflection,  5,  pitch  rate,  q,  lead  angle,  X,  the  relative 

line  of  sight  of  the  target ,  E,,,  ,  tracking  error  e ,  attacker 

A 

normal  acceleration,  A  ,  and  target  motion,  as  well  as  the 

nA 

squares  of  all  of  these  variables.  An  example  of  a  typical 
run  as  recorded  on  the  strip  chart  is  shown  in  Figure  12.' 

Data  Recorders 

Since  the  target  motion  is  filtered  Gaussian  white  noise 
and  the  system  is  linear,  the  statistics  of  the  system  states 
are  assumed  to  be  stationary  and  ergotic  and  the  mean  and 
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Figure  12*. 


Typical  Tracking  Performance i 
Range  -  1000  Feet,  5  G  Target. 


TARGET  ATTACKER  TRACKING  ERROR  LINE  OP  SIGHT 

ACCELERATION  (G*  s)  ACCELERATION  (G's)  (DEGREES)  (DEGREES) 


GA/MA/84M-1 


20.0 

10.0 

0 

-10.0 

-20.0 


0  20  40  60  80  104 

TIME  (SECONDS) 


Figure  12b.  Typical  Tracking  Performance! 

Range  -  1000  Feet,  5  G  Target 
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autocorrelation  functions  can  be  approximated  by  (Ref  16:82) 


1  T 

E{x>  *  |  /fl  x  dt  (38) 

l  T 

E{x2}  2  i  /0  x2  dt  (39) 

The  variance  is  then 

ox2  *  E{(x  -  E{x}]2}  (40) 

*  E{x2}  -  [E{x} J 2  (41) 


The  RMS  value,  or  standard  deviation,  of  each  state  is 
then  the  square  root  of  the  variance.  Thus,  after  each 
100-second  run,  the  outputs  of  the  fourteen  aforementioned 
integrators  were  recorded  and  properly  manipulated  using 
equations  (38) ,  (39) ,  and  (41)  to  determine  the  standard 
deviation  of  each  state  measured.  These  were  recorded  for 
later  comparison  with  the  digital  pilot  analysis. 

Figure  13  shows  the  analog  circuitry  used  to  evaluate 
equations  (38)  and  (39) .  The  potentiometers  were  needed 
for  scaling  purposes  to  prevent  amplifier  overloads. 

The  complete,  tabulated  results  of  the  analog  simulation 
are  attached  to  this  report  as  Appendix  B.  A  discussion  of 
the  results  is  contained  in  Section  V.  A  summary  of  the 
flying  experience  of  the  pilots  who  flew  the  simulation  as 
well  as  some  of  their  comments  pertinent  to  this  study  is 
included  as  Appendix  C. 
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IV.  Analytic  Solution  Using 
an  Optimal  Model  of  the  Human  Controller 

In  the  past,  attempts  to  model  human  controller  be¬ 
havior  have  centered  around  quasi-linear  describing  function 
models  consisting  of  a  lag-lead  network  with  remnant  (Ref  14) . 
Though  these  efforts  were,  for  the  most  part,  successful 
where  the  task  involved  compensatory  tracking,  adaptations 
of  the  model  to  pursuit  tracking  tasks  like  the  one  presented 
here  were  not  as  well  received  (Ref  15).  Work  by  Kleinroan, 
Baron  and  Levison  has  led  to  the  development  of  a  human  con¬ 
troller  model  based  upon  optimal  control  theory  and  the 
assumption  that  a  "well-trained  human  operator  behaves  in 
an  optimal  manner  subject  to  his  inherent  limitations  and 
constraints"  (Ref  2:5). 

The  Kleinman  Model 

An  Overview.  Figure  14  is  a  functional  description 
of  the  optimal  human  operator  model  and  its  interface  with 
system  dynamics.  It  is  assumed  that  the  various  sources 
of  human  randomness,  or  remnant,  are  manifested  as  errors 
in  perception  of  displayed  variables  and  in  execution  of 
desired  control  movements.  Since  the  model  is  linear, 
these  errors  are  lumped  as  observation  noises,  v^,  and  motor 
noise,  vu,  respectively.  These  and  vu  noises  are 
assumed  independent,  Gaussian  noises  of  sufficient  band¬ 
width  as  to  be  considered  white-noise  processes  with 
autocovariances 
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Figure  14,  Functional  Schematic  of  the  Optimal  Pilot  Model  in  the  Control  Loop. 
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B<Vt>VyTU>>  =  V<t_T> 

(42) 

E{vu(t,vu(T,}  =  V(t"T) 

(43) 

The  inherent  time  delays  characteristic  of  human 
controller  behavior,  such  as  perception  and  reaction 
delays,  are  also  lumped  as  an  equivalent  perceptual  time 
delay,  t.  Typical  values  of  t  are  0.2  ±  0.05  seconds 
(Ref  13; 30) .  This  delayed,  noisy  observation  is  operated 
on  by  a  Kalman  estimator  in  cascade  with  a  least  mean  square 

A 

predictor  to  yield  a  "best  estimate",  x(t)  of  the  state  at 
time  t,  conditioned  on  the  observed  output  £(T),T<t.  This 
estimate  is  then  weighted  by  a  set  of  optimal  feedback  gains, 
-Jt*,  to  produce  a  commanded  control  input,  u  (t)  .  Note 

MM  W 

that,  for  this  work,  the  control  input  is  scalar. 

Limitations  on  the  pilot's  ability  (or  willingness)  to 
produce  rapid  or  excessive  control  movements  is  accounted 
for  ny  a  first-order  lag  of  the  form  (t^s  +  l)”1.  In  the 
model,  this  lag  acts  upon  the  noisy  commanded  control  input 
to  produce  the  actual  control  input  to  the  system  dynamics, 
u(t) .  This  lag  effect  creates  the  need  for  defining  u(t) 
as  the  (n+l)th  state  variable.  This  gives 


t^u  ( t)  +  u  ( t )  =  u(t)  =  uc(t)  +  vu(t) 


(44) 


or 


i(t)  5  Vi  -  -  *n+i +  k;  uc(t>  +  k;  V« 


(45) 


Where  u  (t)  is  now  defined  as  the  control  input  to  the 

j 

augmented  system  dynamics. 


35 


GA/MA/74M-1 


Creating  the  Optimal  Control .  The  original  system 
model,  assumed  linear  and  time-invariant,  is  given  by 


X(t)  =  AX (t)  +  b  u(t)  +  w (t) 

(46) 

*(t)  -  H  x  (t) 

(47) 

where  £<t)  is  the  vector  of  observed  variables  and  w(t) 
is  a  zero-mean,  Gaussian,  white  noise  vector  with  co- 
variance 

W  =  E{w(t)wT(T)>  =  W  5 ( t— t )  (48) 

The  augmented  system  is  now  described  by 


z(t) 

=  Aj^eft)  -  b^  uc(t)  +  Wj^t) 

(49) 

z(t) 

«=  col  [x  (t)  ,  u ( t)  ] 

(50) 

Wj^t) 

=  col  [w( t)  ,  vu(t)/xN) 

(51) 

II 

H 

Xi\ 

CO  1(0,  •  •  •  ,0,  1/ Tjj] 

(52) 

and 


(53) 


(54) 


Ai  - 


A  I  b 
_  J _ 

0  |  -l/xN 


The  covariance  of  w^(t)  is  now  given  by 


Efw^tjw^  (t)  }  = 


W  '  0 

I 


0  'VV 
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The  function  of  the  pilot  model  is  to  determine  the 
control  input  u(t)  that  minimizes  the  quadratic  cost 
functional 


J(u)  = 


lim 

r+« 


T 

E {/  (xTQx  +  ru2 
0 


+  gu2)  dt} 


n 

l 


n 

l 


i=l  j=l 


lij  Etx.Xj)  +  rou2  +  go^2 


(55) 


conditioned  on  the  observation  jr(t) .  In  other  words,  the 
model  is  to  minimize  a  weighted  combination  of  the  co¬ 
variances  of  selected  system  states  and  the  mean  squares  of 
control  input  and  control  rate.  It  can  be  shown  that  the 
optimal  control,  u(t) ,  is  generated  by  the  linear  feedback 
law  (Ref  12:360-362). 


tNu(t)  +  u(t)  =  uc(t)  +  vu(t)  =  -t*x  +  vn(t)  (56) 


where 


tN  1/An+1 


l*i  =  TNAi'  1  * 


u  (t)  =  -f*x 

C  —  — 


(57) 

(58) 

(59) 


and  the  optimal  gains,  l*,  and  the  time  constant,  t^,  are. 

identical  to  those  that  are  obtained  from  the  solution  of 

the  noise-free  regulator  problem,  i.e.,  vu=0  (Ref  2:13). 

The  n+1  feedback  gains  are  computed  from  (Ref  12:360-362). 

T 

l  =  (60) 
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Where  Kq  is  the  unique,  positive  definite  solution  to  the 
n+1  dimensional  Riccati  equation 


A  *K  +  K  A  +  Q  -  K  bb  TK  /g  *  0 
o  o  o  o  o  o — o—o  o 

A  and  b  are  the  augmented  matrices 
o  — o 


“1 

> 

It T 

"  0  “ 

i 

• 

- 1 - 

;  So  * 

• 

0  1  0 

6 

1 

1 

and  Qq  is  the  augmented  matrix 


Q 

1 

0 

i 

— 

-1- 

-  — 

0 

1 

0 

1 

(61) 


(62) 


The  Inner  Workings .  The  pilot  is  presented  a  linear 
combination  of  the  state  by  means  of  a  display,  the  output 
of  which  is  modelled  by 

Z'(t)  =  i^zft)  (63) 

where  is  the  constant,  augmented  observation  matrix 
HjMHsOl,  and  (t)  =col[^(t)  :0] .  It  is  assumed  that  the 
human  operator  can  perceive  not  only  the  state  information 
displayed,  but  also  (and  only)  the  first  derivative  of  that 
information  (Ref  12:359) .  Due  to  such  remnant  effects  as 
random  perturbations  in  human  response  characteristics, 
random  errors  in  interpretation  of  the  display,  and  noisy 
display  systems,  the  pilot  actually  perceives  a  delayed. 
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noisy  rendition  of  the  true  system  output, 

^p(t-T)  =  HjZ/t-T)  +  Vy(t-T)  (64) 

where  t,  the  equivalent  perceptual  time  delay  of  the  model, 
defined  earlier,  was  chosen  as  0.2  seconds  for  this  study. 
From  this  perceived  output,  the  Kalman  filter  produces 

A 

a  least  mean  squared  estimate  z(t-t)  of  the  delayed  state 
z(t-x)  through  the  solution  to  the  differential  equation 
(Ref  2:12) 

z(t-T)  *  (Ax  -  CH1TVy-1  H1]  z(t-T) 

+  CHXTV  (t-T)  +  b1uc(t-x)  (65) 

where  C,  the  error  covariance  matrix  satisfies,  in  the  steady 
state  (Ref  12:362). 

AjC  +  CAj1  +  Wx  -  CH^Vy"1  HjC  =  0  (66) 

The  predictor  then  generates  a  current-time  estimate 

A  A 

£(t)=col[x,u(t) J  from 

A  T  ^ 

z(t)  =  e  1  Iz(t-T)  -  ;(t-T)]  +  £(t)  (67) 

where 

£(t)  =  A^(t)  +  bx  uc(t)  (68) 

Figure  15  is  a  detailed  diagram  of  the  computational 
flow  involved  with  this  model.  Recall  that  the  motor  noise 
present  in  Figure  14  is  now  combined  with  the  noise  input 
w^(t)  to  the  system  dynamics.  Also  note  that  in  the 
diagram  as  defined  at  the  n+1  dimensional  row  vector  [£*:0J. 
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Figure  15.  Computational  Flow  Diagram  of  the  Optimal  Pilot  Modal  in  the  Control  Loop  (System 
Equations  Augmented). 
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Finally,  the  system  states  statistics  are  generated 
in  the  form  of  the  covariance  of  £(t)  through  the  expres¬ 


sion  (Ref  12:333) 


T  T 

Mt  ii  T  A.  x  x  A.  0  A.  o 

2  *  E{z(t)zA(t) }  -  e  A  Ce  A  +  /  e  x  W.e  x  do 
—  —  OX 


+  /  eAoe^lTCH. TV  ” * H. Ce^1  TeAo  do  (69) 
o  x  y  i 


where 


I  - 

-*^TN  <  “1/fTN 


The  solution  to  this  equation  can  be  obtained  in  closed- form 
by  noting  that  expressions  of  the  form 

00  — .  — 

X  *  /  e  At  Y  eA  t  dt  (71) 

o 

can  be  evaluated  from  the  solution  to  (Ref  2:143) 

AX  *  XA1  +  Y  *  0  (72) 

Further  computational  methods  for  the  evaluation  of  (69)  are 
described  in  Reference  11.  The  standard  deviations,  or  RMS 
values,  of  the  system  states  are  then  computed  as 


Application  to  the  Air-to-Air  Task 

The  State  Equations  of  Motion.  In  the  introduction 
to  Section  III  the  dynamic  variables  describing  the  air-to- 
air  task  were  listed.  They  were  repeated  here  as  they  were 
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defined  for  the  state  variable  format 


Also, 


X,  =  dummy  state 


X8  '  TT 
*9  5  9 


Vs 


Equations  (10),  (11),  (20),  (27),  (29),  (31),  and  (35) 
together  with  the  noise  filter  equations  are  represented 
in  the  system  dynamics 

x(t)  =  A  x(t)  +  b  u(t)  +  w(t) 


(74) 


in  which  the  "A"  matrix  is 
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Also,  b  =  col  IK_/t  00000000] 

■"  b  a 

and  w(t)  =  col  [0  0  0  0  0  £  0  0  0  ] 

where 


4;  =  white,  Gaussian  driving  noise 

+  This  term  is  zero  for  all  cases  studied  in  this  report. 
*  dummy  feedback  on  this  state  for  stability  precautions. 
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Pros  the  display  the  pilot  can  extract  the  computed 
lead  angle,  X,  the  relative  line  of  sight  to  the  target, 

E_  ,  and  error,  e,  plus  the  rates  of  change  of  each  of  these 

A 

quantities.  However,  E_  ,  X,  and  e  are  a  linearly  dependent 

aA 

set,  so  that  only  two  of  these  variables  need  be  represented 
in  the  observation  vector,  £(t) .  The  error  and  line  of 
sight  were  used  along  with  their  respective  rates  because 
it  was  judged  that  the  pilot  was  more  consciously  aware  of 
these  quantities  than  the  lead  angle.  Thus  the  display  is 
modelled  by 

£<t)  =  Hx (t)  (75) 

with 


and  H  being 
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Cost  Functional  Weightings.  The  weightings  on  the 
terms  of  the  cost  functional,  equation  (55),  were  chosen 
with  the  task  in  mind  and  after  discussions  with  the  pilots 
who  flew  the  simulation.  First,  since  only  fighter  aircraft 
would  be  involved  in  this  task  and  failure  of  the  pilot  to 
perform  well  could  have  catastrophic  results,  the  weighting 
on  control  input,  r,  was  taken  to  be  zero.  In  this  situation, 
the  pilot's  control  inputs  would  only  be  limited  by  the  MG" 
loadings  he  could  physically  withstand.  Of  course,  theoret¬ 
ically  he  could  withstand  an  infinite  number  of  *G's"  in 
this  fixed-base  simulation. 

When  performing  air-to-air  tracking,  the  pilot  is  aware 
primarily  of  the  tracking  error  and  the  rate  of  change  of  that 
error.  Both  are  kept  to  a  minimum  by  the  pilot;  the  former 
for  obvious  reasons,  and  the  latter  because  excessive  rates 
lead  to  overshoots  and  pilot-induced  oscillations  (PIO)  that 
lead  to  even  greater  error.  Further,  the  emphasis  on  low 
tracking  error  rate  is  a  consequence  of  the  "smoothness"  of 
his  flying  technique  that  is  inbred  in  a  pilot  from  his  first 
flying  lesson.  Therefore,  it  was  decided  to  weight  the  error 
and  its  rate  equally  in  the  cost  functional.  Since  error  is 
not  a  system  state,  but  is  a  linear  combination  of  states  [see 
equation  (37)]  in  the  system  output  equation,  a  weighting  matrix 


Q'  in  X  Q  Z  where 


Q"  = 


1 

0 

0 

0 


0 

1 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 


(76) 
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T  T 

would  yield  the  same  effect  as  Q  in  x  Qx  when  Q  =  H  Q'H. 

Finally,  the  weighting  on  control  rate,  g,  is  specified 
when  the  neuro-motor  lag  time  constant,  t^,  is  specified. 

This  is  due  to  the  relationship  between  tn  and  1R+^  given 
by  equation  (57)  and  the  fact  that  the  solution  to  the  Riccati 
equation  (61)  relies  upon  a  knowledge  of  g.  Emperical  data 
have  shown  that  xN  is  of  the  order  of  0.1  to  0.6  seconds, 
with  =  0.1  seconds  being  typical  (Ref  12:363).  Conse¬ 
quently  the  0.1  value  was  used  in  this  study. 

Determination  of  the  Noise  Covariances.  Since  the 
observation  and  motor  noises  depend  upon  the  quality  of  the 
display,  distractions  related  to  the  environment,  and  human 
randomness,  determination  of  numerical  values  for  their 
covariances,  Vy,  and  Vu,  respectively,  can  be  quite  difficult. 
However,  it  has  been  found  that  each  white  observation  noise 

v  (t)  has  a  covariance  that  is  related  to  the  variance  of  its 
Yi 

associated  variable  y^  by  (Ref  12:363) 

Vy  -  irpiE{yi2}  (77) 

where,  on  the  average,  p^=0.01,  which  was  used  herein. 

Similarly,  the  motor  noise  covariance  is  assumed  to 
be  related  to  the  variance  of  the  commanded  control,  u (t), 
by  (Ref  12:363) 

V  =  irp  E{u  2 }  (78) 

u  u  c 

Data-matching  methods  described  in  Reference  12  yield  a 
typical  valve  for  py  of  0.003  which  was  also  utilized  in 
this  study. 
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The  state  equations  for  the  target  motion  noise  filters 


where 


x6  =  ”  tt  x6  +  K 


x?  =  x6  -  x? 


x"  =  Ax"  +  w"(t) 


i  X  ■ 


w"(t)  » 


w6(t>- 


w?(t) 


Note  that  the  "primes"  indicate  subma trices,  not  the  transpose 


The  covariance  of  w"(t)  is 


W'  =  E{w"  (t)  w"A  (t) }  • 


E{£2} 


The  evaluation  of  W"  is  made  from  (Ref  11:1-6) 
k'X'  +  x'a't  +  W"  =  0 


where  X"  is  the  symmetric  covariarce  matrix  of  x"(t) 
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X  =  E{x'(t)  x'T(t) } 

X66  X67 

X76  X77 

and 

Xgg  *  0  x6  =  variance  of  xg(t) 
X77  =  0  x7  =  variance  of  x7(t) 
X?6  »  Xg7  =  E{x6(t)x?(t)} 


(86) 


Since  the  variance  of  x7  is  just  the  square  of  the  RMS  target 
accelerations  stipulated  in  the  simulation,  equation  (85) 
can  now  be  evaluated  to  obtain 


o 

o 


(87) 


where 


Now,  Xg  is  the  only  state  affected  directly  by  w(t) ,  so  W, 
the  covariance  of  w(t)  is  null  except  for  element  wc_. 


Modelling  Threshold  Effects 

In  any  visual  task  there  exists  a  limit  on  the  minimum 
arc  of  resolution  that  the  controller  can  physically  perceive. 
This  is  known  as  the  visual  threshold.  Also,  there  usually 
exists  a  minimum  error  signal  to  which  a  human  may  not  respond. 
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This  is  known  as  the  indifference  threshold.  The  former 
threshold  can  be  qualitatively  evaluated  through  data-matching; 
however  the  latter  would  have  a  high  dependency  on  the  situation 
involved.  For  instance,  a  pilot  may  not  be  too  concerned  with 
a  certain  glide  path  deviation  at  five  miles  from  touchdown 
in  the  landing  approach,  whereas  that  same  deviation  may  be 
totally  unacceptable  at  one  mile  range. 

In  the  pilot  model  for  this  task  these  effects  are  not 
distinguished  and  are  modelled  as  a  statistical  linearization 
of  the  "deadzone"  nonlinearity  shown  in  Figure  16  where 


(Ref  9:87-95) 


f  (y) 


P 

y  -  a,  y  >  a 
0  ,  -a  <  y  <  a 

y  +  a,  y  <  -a 


(88) 


This  amounts  to  replacing  this  nonlinearity  with  an  equivalent 

A  • 

gain  f(y).  This  gain,  always  less  than  unity,  is  plotted 
versus  a */ (ST~  a x)  in  Figure  16,  as  taken  from  Ref.  6,  page  238. 
Here  the  parameter  "a"  is  the  actual  minimum  deviation  that 
can  be  resolved  by  the  observer.  It  has  been  found  that  0.05° 
is  the  typical  minimum  arc  of  resolution  (Ref  9:93).  Since 
in  this  simulation  20  inches  separated  the  observer  from 
the  display,  this  0.05°  converts  to  approximately  0.05cm 
on  the  oscilloscope.  To  accomodate  the  large  lead  angles 
generated  at  the  3000- foot  range,  it  was  necessary  to  scale 
the  display  to  approximately  13°  per  cm.  Therefore,  the 
error  threshold,  a£,  for  this  simulation  was  computed  to 
be  0.65°.  Typical  values  for  rate  threshold  are  ~0.05  to 
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Figure  16.  Threshold  Nonlinearity  and  Equivalent  Gain  Plot 
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0.10  degrees  of  visual  arc  per  second  (Ref  8:15).  Using 
the  upper  limit  in  this  simulation,  the  error  rate  threshold, 
a*,  was  computed  to  be  1.3  degrees  per  second. 

C* 

Thus  the  image  presented  to  the  pilot  appears  as 

y  (t)  *  ^(yJy^t)  +  v  (t)  (89) 

pi  11  yi 

or 

y*  (t)  -  y. (t)  +  fT‘(y)  v  (t)  (90) 

PA  i  1  yi 

Therefore,  when  the  effects  of  threshold  are  taken 
into  account,  the  covariance  of  each  observation  noise  is 
given  by 

V  =  up  fT2(y)  E{y?}  (91) 

yi  1  1  i 

The  Pilot  Model  in  Action 

The  objective  of  this  chapter  has  been  to  present  the 
equations  and  rationale  of  the  Kleinman  pilot  model  applied 
to  the  air-to-air  tracking  task.  In  practice,  the  compu¬ 
tations  are  performed  by  means  of  a  digital  computer  program. 
The  program  used  in  this  study  was  adapted  from  that  pre¬ 
sented  in  Reference  13  by  Major  James  D.  Dillow,  Assistant 
Professor  of  Mathematics  at  the  Air  Force  Institute  of 
Technology,  Wright-Patterson  Air  Force  Base,  Ohio.  A  brief 
description  of  that  program  follows. 

After  the  situational  data,  e.  g.,  stability  derivatives, 
velocities,  altitudes,  target  noise  covariance,  etc.,  and  cost 
functional  weightings  are  input  to  the  program,  the  feedback 
gains,  l,  and  consequently  tn,  are  computed.  If  the  desired 

value  of  x.,  is  not  achieved,  the  weighting  on  the  control 
N 
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rate,  g,  is  adjusted  and  the  computations  repeated.  This 
process  is  repeated  until  the  proper  weighting  has  been 
discovered  that  yields  the  desired  t^.  The  optimal  feedback 
gains  are  now  known. 

The  uncorrelated  observation  and  motor  noises  are  now 
added  to  the  system  in  the  form  of  initial  guesses  at  the 
values  of  their  respective  covariance  matrices  Vy  and  Vy. 

Enough  information  is  now  known  to  compute  the  error  covariance, 
C,  from  equation  (66).  The  covariances  of  the  system  states, 
which,  in  the  augmented  form,  include  the  control,  u(t) ,  are 
then  computed  from  equation  (69) .  Finally,  new  values  for 
Vy  and  Vu  are  calculated  using  equations  (91)  and  (78) , 
respectively.  If  the  differences  between  the  new  and  old 
values  of  the  noise  covariances  are  not  tolerable,  the  process 
is  repeated  until  convergence  on  the  correct,  within  toler¬ 
ances,  covariances  is  achieved.  The  RMS  values  of  the  pilot 
model  performance  indexes  are  then  computed  from  equation  (73) . 

The  results  obtained  from  the  pilot  model  "flying" 
all  twelve  test  cases  is  compared  in  the  next  section  with 
those  obtained  from  the  analog  simulation. 
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V.  Analysis  of  Results 

Table  II  and  Figures  17  through  21  are  a  summary  of 
the  preliminary  findings  of  this  study.  There  is  very 
definitely  a  one-to-one  correlation  between  predicted  and 
actual  elevator  deflection,  lead  angle,  and  line  of  sight. 
However,  it  would  appear  at  first  glance  that  the  pitch 
rate  and  error  plots  do  not  correlate  very  well — a  very 
serious  flaw,  since  the  ability  to  track  with  minimum  error 
is  a  prime  consideration  when  judging  the  effectiveness  of 
a  weapon  system  in  this  task.  A  second  look  is  more  re¬ 
vealing.  Look  first  at  pitch  rate.  In  every  case,  the  com¬ 
puter  generated  more  pitch  rate  than  did  the  pilots  in  the 
simulation.  However,  there  is  a  marked  distinction  between 
pitch  rate  magnitudes  at  the  3000 -foot  range  and  those  at 
1000  feet.  For  clarity,  a  line  has  been  drawn  separating 
the  two  groups  on  the  plot. 

Now  look  at  the  plot  of  tracking  eivor .  Though  the 
data  points  are  somewhat  scattered,  in  the  mean  there  does 
appear  to  be  almost  a  one-to-one  correlation  between  pre¬ 
dicted  and  actual  error  at  the  1000 “foot  range,  the 
correlation  becoming  better  as  the  target  motion  reduces 
from  5.0G  to  3.5G.  The  3000-foot  data  is  another  story. 

A  line  drawn  through  the  mean  of  this  grouping  has  a  slope 
of  almost  one-half,  meaning  that  the  computer  is  tracking 
nearly  twice  as  well  as  the  simulation  pilots. 

The  pitch  rate  and  error  grouping  tendency  led  to  a 
deeper  analysis  of  the  differences  in  dynamics  at  the  two 
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Figure  17.  Predicted  Versus  Actual  RMS  Elevator  Deflection  (Degrees) 
Equal  Weighting  on  Error  and  Error  Rate. 
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Figure  21.  Predicted  Versus  Actual  IMS  Error  (Degrees) 
Equal  Weighting  on  Error  and  Error  Rate. 


ranges.  The  answer  was  found  in  the  sight  equation,  (20). 

As  was  stated  in  Section  II,  in  practice  the  sight  parameters 
and  are  normally  held  constant  at  some  nominal  value 
determined  from  nominal  conditions  of  altitude,  temperature, 
velocities,  and  range.  However,  in  this  study  the  sight 
time  constant,  T^,  the  projectile  time  of  flight,  and  V^, 
the  projectile  average  relative  velocity,  were  recomputed 
for  each  case  by  the  method  shown  in  Appendix  A.  In  effect, 
then,  changing  ranges  changed  the  sight,  just  as  if  the 
"black  boxes”  in  the  aircraft  had  been  changed!  The  difference 
is  amazing.  At  1000  feet,  T^  averaged  0.34  seconds.  This 
compares  with  1.38  seconds  at  3000  feet--a  factor  of  four. 

Thus,  in  addition  to  the  generation  of  larger  lead  angles, 
the  sight  dynamics  were  accentuated  in  the  form  of  more 
visible  oscillations.  The  pilots  found  that  they  could 
not  be  as  aggressive  at  this  range  because  of  the  sight 
sensitivity.  Large  pitch  rates  generated  large  sight  over¬ 
shoots,  and  the  resulting  long  settling  time  deprived  the 
pilot  of  the  exact  knowledge  of  his  tracking  error.  So 
the  pilots  compensated.  They  reduced  their  pitch  rate  and 
concentrated  on  "easing"  the  sight  onto  the  target.  This, 
of  course,  led  to  longer  periods  in  which  the  sight  was  not 
on  the  target  and,  consequently,  greater  RMS  tracking  error 
than  would  have  been  expected. 

So,  through  subjective  questioning  of  other  fighter 
pilots  it  was  decided  that  relative  weightings  of  e  =  3, 
e  =  1,  and  X  =  2  might  work.  Since  the  cost  functional  is 
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quadratic,  Q'  becomes 


Q'  = 


9 

0 

0 


0 

5 

0 


0 

0 

0 


0 

4 

0 


0  4 


0 


4 


(92) 


With  these  weightings,  the  results  tabulated  in  Table  III 
were  obtained.  Plots  of  predicted  versus  actual  parameters 
are  presented  in  Figures  21  through  26.  Only  slight  improve¬ 
ment  can  be  seen  in  6,  A,  and  I_  .  However,  marked  improve- 

A 

ment  can  be  seen  in  both  pitch  rate  and  tracking  error.  The 
addec  xghting  on  sight  rate  and  the  reduced  weighting  on 
error  rate  has  had  the  effect  of  considerably  increasing 
the  tracking  error  at  3000  feet  while  having  a  smaller  det¬ 
rimental  effect  on  the  1000- foot  errors.  The  data  points  are 
still  far  from  a  perfect  match,  but  they  are  quite  acceptable. 
More  work  needs  to  be  done  toward  finding  the  rationale  for 
computing  an  optimum  set  of  weightings. 

The  results  not  only  correlate  quite  well,  but  they 
also  make  good  intuitive  sense.  Note  that  for  a  given 
range  and  target  motion  the  data  point  out  the  relative 
strengths  of  the  different  aircraft  in  performing  this 
task.  For  instance,  at  the  1000- foot  range  the  F-5  was 
predicted  to  be  the  best  aircraft  for  the  task  (per  "G" 
level),  followed  closely  in  ability  by  the  F-4E,  and  then 
by  the  A-7.  With  one  slight  anomaly  at  the  5G  level,  the 
simulation  data  indicates  likewise.  Since  the  F-5  is 
primarily  an  air-to-air  fighter,  the  F-4E  mere  irulti- 
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purpose  in  design,  and  the  A-7  primarily  an  air-to-ground 
weapon,  the  results  make  sense.  On  the  other  hand,  at  the 
3000- foot  range  the  F-?E  is  judged  the  best,  followed  by 
the  A-7.  Last  is  the  F-5.  The  principal  reason  for  this 
is  the  extremely  responsive  nature  of  the  longitudinal 
dynamics  of  the  F-5.  It  if,  so  pitch  sensitive  that  it 
aggrevates  the  problems  with  the  sight  dynamics  at  this 
long  range.  On  the  other  hand,  the  A-7  is  very  sluggish 
(comparatively  speaking)  in  the  longitudinal  axis,  a  feature 
that  was  very  helpful  in  damping  the  sight  oscillations  at 
this  range. 

It  should  be  noted  that  the  relative  performance 
figures  given  here  are  not  strictly  indicative  of  one  air¬ 
craft's  ability  to  perform  relative  to  another.  The 
situational  parameters  chosen  for  this  work  were  done  with 
an  attempt  in  mind  toward  obtaining  a  good  "spread"  in  the 
data.  The  F-5  at  *5.135  feet  and  0.81  Mach  represents  a 
good,  responsive  environment;  the  A-7  at  15,000  feet  and 
0.61  Mach  represents  a  sluggish  set  of  parameters,  by 
comparison;  and  finally,  the  F-4E  at  15,000  feet  and  0.9 
Mach  would  be  classified  as  a  good  median  situation.  This 
analysis  is  born  out  by  the  groupings  of  data  points  on 
elevator  deflection  in  Figures  17  and  22.  When  these 
groupings  are  compared  with  the  resulting  pitch  rate 
statistics  in  Figures  18  and  23,  it  can  be  seen  that  the 
F-4E  required  somewhat  more  elevator  travel  than  the  F-5 
to  produce  roughly  the  same  amount  of  pitch  rate.  The 
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difference  is  even  more  pronounced  with  the  A-7 .  The 
higher  pitch  rates  generated  with  the  A-7  are  a  consequence 

t 

of  the  response  lag  produced  by  its  sluggish  environment. 

By  the  time  the  aircraft  responded,  tracking  error  was 
greater  and  higher  rates  were  required  to  reduce  that 
error . 

Another  look  at  the  elevator  deflection  plots  indicates 
that,  regardless  of  the  level  of  performance  achieved  in 
any  one  situation,  the  pilots  certainly  had  to  work  much 
harder  with  the  A-7  than  with  the  F-5  or  F-4E.  This  con¬ 
clusion  was  substantiated  by  the  pilots  themselves. 

A  final  note  before  concluding:  there  should  be  no 
attempt  made  to  correlate  aircraft  tracking  performance 
differences  between  the  two  ranges  used.  Since  sight 
dynamics  were  changed  each  time  range  was  changed,  one 
would  be  comparing  an  aircraft  having  a  sight  designed 
for  3000  feet  and  being  used  at  3000  feet  against  one  having 
a  sight  designed  for  1000  feet  and  being  used  at  that  range. 

Obviously,  with  given,  constant  sight  parameters,  tracking 
error  would  almost  always  be  less  at  the  3000- foot  range 
than  at  1000  feet. 
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VI.  Concluding  Remarks  and  Recommendations 
Conclusions 

The  "scatter  diagrams”  of  Section  V  speak  for  them¬ 
selves.  In  fact,  when  applied  to  the  results  of  this  study, 
the  term  "scatter"  is  almost  a  misnomer.  One  can  but  conclude 
that,  given  the  ideal  simulation  conditions  used  in  this 
study  (two-dimensional,  linear,  stationary,  Gaussian,  zero- 
mean)  ,  the  model  of  the  human  controller  exhibited  herein 
is  a  superb  imitation  of  a  human  pilot  when  performing  the 
air-to-air  tracking  task  with  an  LCOS  system.  It  would  appear 
that  one  need  only  give  judicious  thought  toward  evaluation 
of  the  proper  weightings,  Q',  on  the  displayed  variables  and 
toward  consideration  of  the  threshold  effects  of  the  display 
available  in  order  to  be  able  to  predict  the  performance 
of  given  aircraft,  flight  control,  and  sight  combinations. 

The  only  weak  link  in  the  model  seems  to  lie  in  the 
selection  of  these  weightings  on  state.  They  are,  of  course, 
critical  to  the  success  of  the  predictions  derived  from  the 
execution  of  the  pilot  model.  The  predicted  magnitudes  of 
lead  angle,  elevator  deflection,  and  line  of  sight  were 
found  to  be  fairly  insensitive  to  changes  in  the  weightings. 
However,  pitch  rate  and  error  were  very  sensitive  to  these 
changes  and,  since  tracking  error  is  the  prime  factor  in 
the  evaluation  of  system  performance,  great  care  must  be 
taken  in  the  selection  of  these  weightings.  They  must 
therefore  reflect  intimate  knowledge  of  human  controller 
behavior  in  the  environment  created  for  him  by  the  task 
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at  hand. 

An  alternative  to  rigorous  selection  of  relative  state 
weightings  would  be  to  use  such  weightings  as  are  given  in 
this  report  and  use  the  pilot  model  to  obtain  performance 
data  on  a  system  of  known  performance  quality.  Then  other 
systems  could  be  analyzed  based  upon  performance  relative  to 
this  base  system. 

Recommendations 

The  success  of  this  study  should  open  the  door  to 

expanded  study  of  the  application  of  this  pilot  model  to 

other  tasks,  such  as  air-to-ground  weapon  delivery.  For 

the  present,  however,  there  is  much  continuation  work 

needed  in  the  application  of  the  model  to  the  air-to-air 

task.  The  next  step  should  be  to  perform  a  time-varying 

simulation  which  would  account  for  closure  of  the  attacker 

on  the  target.  This  would  require  extension  of  the  pilot 

« 

model  to  the  time-varying  case.  The  results  of  this  study 
should  then  be  compared  with  the  results  published  herein 
to  determine  if  the  simpler,  stationary  model  is  not 
adequate  for  evaluation  of  system  performance. 

Of  course,  three-dimensional  modelling  and  testing 
of  the  pilot  model  in  the  air-to-air  environment  is  definitely 
warranted.  No  study  is  complete  until  this  has  been  done. 

The  lateral  and  yaw  dynamics  thus  added  would  require 
formulation  of  a  pilot  model  capable  of  multiple  control 
outputs . 

Further  extension  of  the  model  to  include  non linearities 
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such  as  varying  velocities  and  velocity  jump  should  then 
follow.  Once  again,  the  pilot  model  would  have  to  be  modified 
to  accept  a  nonlinear  plant. 

Another  possibility  would  be  the  development  of  a 
model  for  the  error  propagation  for  the  bullet  as  a  function 
of  the  attacker  states  and  target  motion  so  as  to  compute 
the  probability  of  a  kill. 

Finally,  contingent  upon  the  outcome  of  this  later 
research,  industry  should  seriously  consider  the  potential 
of  this  model  for  providing  inexpensive,  much-needed  testing 
of  integrated  subsystems  at  the  "idea"  level.  This  testing 
would  reveal  any  serious  deficiencies  in  required  performance 
characteristics  before  the  actual  hardware  is  built.  In 
fact,  this  approach  could  be  used  to  synthesize  an  "optimal" 
integrated  fire  control/flight  control  system  in  a  relatively 
inexpensive  manner. 
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Appendix  A 


Determination  of  Projectile  Time  of  Flight 


And  Average  Relative  Velocity 


Equation  (20) ,  used  to  compute  the  required  lead 
angle,  X,  is  dependent  upon  two  projectile  trajectory 
parameters,  the  values  of  which  are  not  altogether 
readily  apparent.  They  are  T^,  the  time  of  flight  of 
the  projectile  to  impact  with  the  target,  and  V^,  the 
average  velocity  of  the  projectile  relative  to  the  weapon 
station  over  this  time  of  flight.  Both  parameters  are. 
of  course,  dependent  upon  each  other  as  well  as  target 
range,  velocity,  and  maneuvering  tactics.  This  addendum 
relates  the  simple  manner  in  which  "ball  park"  estimates 
of  these  two  parameters  were  computed.  Perfection  was 
not  the  goal,  so,  with  respect  to  some  of  the  assumptions 
made,  realism  was  sacrificed  somewhat  in  favor  of  simplicity. 

For  the  purposes  of  this  study,  nominal  values  for  T^ 
and  were  computed  based  upon  the  tracking  ranges, 
altitudes,  and  airspeeds  to  be  used  in  the  simulation. 

The  target  was  limited  to  constant  3.5  G  and  5.0  G  turns. 

It  was  assumed  that  the  tracking  aircraft  had  achieved 
the  proper  lead  angle  and  that  projectile  impact  on  the 
target  would  occur  after  T^  seconds  of  flight.  Figure  27 
shows  the  geometry  involved. 
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Figure  27.  Projectile  Trajectory  Geometry 


The  attacker  angle  of  attack  was  assumed  zero  so 

that 

l£il  -  IX*  +  S«l  -  +  I2«l  <”> 

where 

Vj  =  initial  projectile  total  velocity 
=  projectile  muzzle  velocity 
=  attacking  aircraft  velocity 

Current  techniques  in  ballistics  modelling  use  the  integral 
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to  calculate  projectile  velocity,  Vp,  at  any  time,  t,  after 
firing.  Here  p  is  the  atmospheric  density  at  the  flight 
altitude,  pQ  is  sea  level  atmospheric  density,  and  is 
a  ballistic  constant  producing  best  correspondence  with 
appropriate  ballistic  tables  (Ref  21:3).  KQ  was  chosen 
as  0.00614  in  this  study  from  an  analysis  of  ballistics 
tables  F-7  and  F-8  on  page  F-16  of  Reference  5. 

The  range  of  the  target  from  the  position  the  attacker 
occupied  at  the  time  of  firing  is  given  by 


St 


\ 

\ 


(95) 


So, 

the 


from  the  geometry  given  in  Figure  27, 
target  can  be  modelled  by 


the  velocity  of 


VT  cos  yt 
VT  sin  yT 


(96) 


As  was  pointed  out  in  Section  III,  yt  can  be  computed 
from  the  solution  to 


Finally,  the  time  rate  of  change  of  projectile 
position,  Rp,  is 


(35) 


(97) 
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So,  equations  (35),  (94),  (96),  and  (97)  were  integrated 

in  a  digital  computer  program  using  a  Runge-Kutta  integration 

routine.  At  each  step,  the  range  of  the  projectile, |R  |, 

P 

and  the  range  of  the  target,  given  by 


|Rj  =  /  RT  2  +  RT  2  (98) 

x  y 

were  compared  to  determine  if  impact  had  occurred.  When 
impact  occurred,  the  elapsed  time  in  the  Runge-Kutta  routine 
was  then  T^.  The  average  velocity  of  the  projectile  was 
then 


1  Tf 

V  =  =r-  f  V  dt 
av9  Tf  o  p 


It  can  be  shown  that  (Ref  7:66) 

V  =  VT  e  p 
P  I 


Therefore, 


VI  Tf  -VRo 
i  =  —L  /  r  e  Bpdt 
'avg  Tf 

VT  — K  pR 

=  _ ?_  f!  -  e  B  P) 

K_pR  1  e  > 

BH  P 


The  average  relative  velocity  of  the  projectile,  V^, 
after  T^  seconds  of  flight  v.’as  then 


(99) 


(100) 


(101) 

(102) 


V 


f 


V. 


(103) 
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Tabulated  Simulation  Data 


Table  IV, 
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Appendix  C 

Pilot  Data  And  "Pireps"* 

The  three  pilots  that  "flew"  the  simulation  were 
chosen  for  their  flying  experience  in  fighter- type  aircraft. 

Table  VIII  summarizes  that  experience.  All  three  are 
presently  pursuing  master  of  science  degrees  in  engineering 
at  the  Air  Force  Institute  of  Technology  and  are  thus  well 
qualified  to  comment  on  the  realism  and  authenticity  of  the 
simulation. 

Given  the  fact  that  only  longitudinal  dynamics  were 
simulated,  the  pilots  generally  thought  that  the  simulation 
quite  adequately  represented  vehicle  and  sight  dynamics. 

In  fact,  it  was  suggested  that  maybe  this  simulator  would 
by  an  inexpensive  means  of  familiarizing  new  fighter  pilots 
with  the  peculiarities  of  LCOS  dynamics,  and  how  to  handle 
them. 

Principal  among  the  complaints  concerning  the 

simulation  was  the  fact  that  extensive  negative  "G's"  were  i 

\ 

necessarily  required  to  track  the  zero-mean  target  motion.  j 

) 

i 

This  just  does  not  happen  in  reality.  However,  after  a  f 

J 

period  of  training,  the  pilots  became  used  to  pushing  on 

the  stick  and  their  performance  was  not  affected.  The  f 

fact  that  there  was  no  way  to  feel  the  aircraft  "G"  loading  f 

i 

in  the  fixed-base  simulator  was  also  deemed  detrimental,  j 

I 

The  term  "pirep"  comes  from  the  weather  forecasting 
terminology  for  pilot  reports  on  weather  conditions 
encountered  aloft. 
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Tab.’e  VIII. 

Stanry  of  Pilot  Experience 


Pilot 

(Total  Plying  Time) 

Type 

Aircraft 

Number  Hours 

This  Aircraft 

Years  Flown 

T-37 

125 

1963  -  1964 

Pilot  #1 

T-33 

325 

1965  £  1971  -  1972 

(2250  Hours) 

F-4C  £  D 

1800 

1965  -  1970 

T-37/T-38 

290 

1966  -  1967 

P-102 

65 

1967 

Pilot  #2 

P-101 

420 

1967  -  1969 

(2100  Hours) 

T-33 

380 

1967  -  1969 

P-4 

400 

1969  -  1970 

Hawker  Hunter 

25 

1970 

B.A.C.  Lightning 

520 

1971  -  1973 

T-37/T-38 

235 

1966  -  1967 

Pilot  #3 

P-100 

680  • 

1967  -  1969 

(1959  Hours) 

F-4 

1044 

1970  -  1973 
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but  not  greatly  important. 

Finally,  the  display  on  the  oscilloscope  did  not 
allow  the  pilot  the  ability  to  see  his  "angle-off" 
the  target — that  angle  formed  by  the  extension  of  the 
longitudinal  axes  of  the  two  aircraft.  Again,  the  degrading 
effect  was  not  deemed  important. 

In  summary,  they  felt  the  simulation  dynamics  were 
quite  authentic,  however  the  negative  "G"  tracking  and 
the  primitive  display  tended  to  detract  from  the  realism 
intended.  It  should  be  noted  that  this  fact  did  not 
dampen  their  desire  to  do  their  very  best  in  the  performance 
of  the  task  presented  them. 
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Vita;  Thongs  R.  Harvey 


Thomas  R.  Harvey  was  born  an  Army  "brat" 


nown  today  as 


Vandenberg  Air  Force  Base.  After  years  of  wandering  from 
fort  to  fort,  his  family  finally  settled  in  Gainesville, 


Florida.  Graduation  froi 


was 


followed  by  four  years  at  the  Virginia  Military  Institute, 
from  which  he  graduated  in  1966  with  a  Bachelor  of  Science 
degree  in  Electrical  Engineering,  a  commission  as  a  Second 
Lieutenant  in  the  United  States  Air  Force,  and  a  ticket  to 
flight  school.  In  1968,  after  successfully  completing 
Undergraduate  Pilot  Training  at  Laredo  Air  Force  Base, 
Texas,  he  upgraded  to  the  F-100  aircraft  and  subsequently 
spent  a  14-month  tour  at  Wheelus  Air  Base,  Libya.  There 
he  was  a  range  safety  officer  on  the  USAFE  air-to-air  and 
air-to-ground  gunnery  ranges,  while  maintaining  his  own 
proficiency  in  the  F-100C  and  U-6A  aircraft.  Upon  his 
return  to  the  United  States  in  1970,  he  was  assigned  to 
fly  0141 A  aircraft  out  of  Dover  Air  Force  Base,  Delaware. 
Two  and  a  half  years  of  flying  global  airlift  ensued  during 
which  time  he  attained  the  position  of  aircraft  commander. 
In  1972,  he  was  selected  for  admission  to  the  Air  Force 
Institute  of  Technology,  which  he  entered  in  August  of  that 
year.  He  graduated  with  the  Master  of  Science  degree  in 
Astronautical  Engineering,  specializing  in  Guidance  and 
Control,  in  March,  1974. 

This  thesis  was  typed  


[Pll  Redacted] 
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